Stickiness in Hamiltonian systems: from sharply divided to hierarchical phase space 
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We investigate the dynamics of chaotic trajectories in simple yet physically important Hamiltonian 
systems with non-hierarchical borders between regular and chaotic regions with positive measures. 
We show that the stickiness to the border of the regular regions in systems with such a sharply 
divided phase space occurs through one-parameter families of marginally unstable periodic orbits 
and is characterized by an exponent 7 = 2 for the asymptotic power-law decay of the distribution 
of recurrence times. Generic perturbations lead to systems with hierarchical phase space, where 
the stickiness is apparently enhanced due to the presence of infinitely many regular islands and 
Cantori. In this case, we show that the distribution of recurrence times can be composed of a sum 
of exponentials or a sum of power-laws, depending on the relative contribution of the primary and 
secondary structures of the hierarchy. Numerical verification of our main results are provided for 
area-preserving maps, mushroom billiards, and the newly defined magnetic mushroom billiards. 
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I. INTRODUCTION 

Hamiltonian systems usually exhibit divided phase 
space, where regular and chaotic regions coexist. An im- 
portant property of chaotic trajectories in divided phase 
spaces is the intermittent behavior with sporadically long 
periods of time spent near the border of regular regions 
Because of this stickiness and the ergodicity of the 
chaotic regions, even small islands can have a large ef- 
fect on global properties of the system, such as trans- 
port and decay of correlations Q . The stickiness can 
be quantified in terms of the distribution P{T) of re- 
currence times T of a typical trajectory to a pre-defined 
recurrence region, usually taken away from regular is- 
lands. For fully chaotic hyperbolic systems, the recur- 
rence time distribution (RTD) decays exponentially 
while for Hamiltonian systems with divided phase space 
the RTD has been argued to decay approximately as a 
power law P{T) ~ T~'^ for large T, where 7' is a scal- 
ing exponent 0, 0, 0, IE 111- For power-law decay, the 
cumulative RTD is given by 



(1) 



T=T 



where 7 = 7' — 1. We say that a system has the prop- 
erty of stickiness if Q{t) decays at least as slowly as t~'^ 
for some 7 > 0. The existence of a finite mean recur- 
rence time implies 7 > 1 0- Experimental evidence of 
stickiness has been observed, for example, in the trans- 
port of particles advected by fluid flows ^ and in the 
fluctuations of conductance in chaotic cavities [a . 
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In Hamiltonian systems, the border between a regu- 
lar and chaotic region often presents a complex hierar- 
chical structure of Kolmogorov-Arnold-Moser (KAM) is- 
lands and Cantori. Cantori are invariant Cantor sets that 
work aspartial barriers to the transport close to KAM 
islands 0, |^. Although many properties of this struc- 
ture are well understood, their consequences to the dy- 
namics are still a matter of intense study Q, |E H ^3 ■ 
Even in the simplest case of two-dimensional systems, a 
number of non-equivalent models have been proposed to 
describe the stickiness of chaotic trajectories. Meiss and 
Ott introduced a Markov-tree model that accounts for 
the hierarchical structure and predicts a scaling exponent 
7 = 1.96 01 . Chirikov and Shepelyansky used renormal- 
ization arguments at the breakdown of the golden mean 
torus to predict a universal exponent 7 = 3 j^. Za- 
slavsky and co-workers applied different renormalization 
arguments to the case of self-similar island chains, ob- 
taining simple relations between 7 and the scaling prop- 
erties of these chains 0. There is also strong evidence 
of other stickiness mechanisms in generic Hamiltonian 
systems 0, 0, 0, . The effects described in these 
previous works typically coexist and are responsible for 
finite-time numerical estimates of 7 lying in the inter- 
val 1.5 < 7 < 2.5 6J. However, because the conver- 
gence in Hamiltonian systems can take an arbitrarily long 
time, in general it is not even clear whether the RTD ap- 
proximates a power-law distribution in the asymptotic 
limit. This slow convergence has inspired Motter and co- 
workers to introduce a new model that accounts for the 
effects of the Cantori structure at finite times [lE| ■ While 
the general asymptotic behavior remains unresolved, the 
insight provided by the study of classes of comprehensible 
Hamiltonian systems is of fundamental importance. 

In this paper, we investigate a mechanism for the stick- 
iness of chaotic trajectories in Hamiltonian systems with 
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non-hierarchical borders between regular and chaotic re- 
gions when both regions can have positive measure. Ex- 
amples of systems with such a sharply divided phase 
space include piecewise-linear area-preserving maps 
ll^ and mushroom billiards [TsL [l9| . In Hamiltonian sys- 
tems, it is a common sense statement to relate the stick- 
iness to the presence of hierarchical structures of KAM 
islands and Cantori. While this statement by itself is 
not wrong, here we show that regular islands with non- 
hierarchical borders also stick. Our results are valid for 
both zero- and positive-measure regular islands, and in- 
clude as particular cases previous findings for systems 
not exhibiting KAM islands, such as the Sinai and sta- 
dium bilhards [lO, HHHIll . We find that the stickiness 
near non-hierarchical borders occurs due to the presence 
of one-parameter families of marginally unstable peri- 
odic orbits (MUPOs). Based on the analysis of MUPOs, 
we show that the recurrence time does follow a power- 
law distribution and that the scaling exponent is 7 = 2 
in two-dimensional sharply divided phase spaces, irre- 
spective of other details of the system. We also study 
the properties of generic perturbations of these systems. 
Based on numerical simulations of mushroom billiards 
perturbed by magnetic fields, we observe that the per- 
turbations generate hierarchical structures of KAM is- 
lands and Cantori of the same nature of those observed 
in typical Hamiltonian systems. The perturbation of 
Hamiltonian systems with sharply divided phase space 
thus represents a new route to Hamiltonian systems with 
hierarchical phase space. Previously considered routes 
start either from fully integrable or fully chaotic config- 
urations. The onset of hierarchical structures introduces 
oscillations in the RTD, which we show to be related to 
the relative contribution of the primary and secondary 
structures of the hierarchy. 

The paper is organized as follows. In Sec. ^ we an- 
alyze the stickiness in sharply divided phase spaces. In 
Sec. liiil we consider the effect of the hierarchical struc- 
tures when a system with sharply divided phase space is 
perturbed. Discussion and conclusions are presented in 
the last section. 



II. SHARPLY DIVIDED PHASE SPACE 

We study the stickiness of chaotic trajectories in two- 
dimensional Hamiltonian systems with non-hierarchical 
borders between the regions of chaotic and regular mo- 
tion. As examples, we consider piecewise-linear area- 
preserving maps and mushroom billiards in Sees. Ill Al and 
IIIBI respectively. In contrast to the previously consid- 
ered stadium and Sinai billiards, in these systems both 
the chaotic and regular regions of the phase space have a 
positive measure. The scaling exponent 7 = 2 is derived 
in Sec. HTn 



A. Piecewise-linear maps 

Consider two-dimensional area-preserving maps of the 
form 

Vn+i = Vn + K f {xn) mod 1, 

Xn+l=Xn+yn+l mod 1, 

where K is & parameter that controls the nonlinearity. 
For f{xn) = sin(27ra;„), Eq. ^ corresponds to the stan- 
dard map, which has served as a prototype of Hamilto- 
nian system in numerous studies of stickiness in hierar- 
chical phase space [H H H Hi] . 

However, for /(a;„) defined as a piecewise-linear func- 
tion of the interval Xn € [0, 1], the phase space of map 
(O can be sharply divided in the sense that regular and 
chaotic regions are separated by a simple curve f^F]. As 
shown in Refs. 2&> 21}, the form and distribution of 
the regular regions in general depend on the function / 
and on the parameter K. In the case of hierarchical dis- 
tribution of islands, it has been shown that the stickiness 
of chaotic trajectories leads to anomalous diffusion in the 
extended phase space of these maps [2^ . 

To quantify the stickiness in the case of sharp border, 
we consider two simple examples with a single regular 
island: 

(i) The first example is obtained for 

/(x„) = l-|2a:„-l|, (3) 

and was called continuous sawtooth map in Ref. 17]. The 
phase space of this map is shown in Fig.^a) for K — 3/2. 
It was argued in Ref. llTll that in this case a single regular 
island exists [Fig. ^a), triangular region]. As we show, 
the absence of other islands and the ergodicity of the 
chaotic region do not rule out the possible existence of 
zero measure sets of MUPOs in the chaotic region. In this 
paper we use the acronym MUPOs to refer to periodic 
orbits in contact with the chaotic component that have 
zero Lyapunov exponents and real eigenvalues with mod- 
ulus one. In sharply divided phase spaces, we regard the 
borders of regular islands as families of MUPOs when- 
ever they are periodic. For example, for the continuous 
sawtooth map, the following sets and their images corre- 
spond to one-parameter families of period-three MUPOs: 
{x = 1/6,1/6 < y < 1/3}, {x = 1/3,1/3 < y < 2/3}, 
and {x — 1/2,1/2 < y < 1}, where the latter is at the 
border of the island. These families of MUPOs corre- 
spond to the straight lines in Fig.^a). 

(ii) A second example of sharply divided phase space is 
obtained for 

C -Xn if < Xn < 1/4, 

fiXn) = { -1/2 + Xn if 1/4 <Xn< 3/4, (4) 
[l-Xn if 3/4 <Xn<l, 

as considered in Ref. |l6j]. The properties of map ©-Q 
with K = 2 are essentially the same of the continuous 
sawtooth map with K = 3/2 [Fig. [Hb) vs. Fig. [Ha)]. 
In particular, the phase space of map ©-lEll has a single 
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FIG. 1: (Color Online) (a) Phase-space portrait of the con- 
tinuous sawtooth map for K — 3/2. The dots correspond 
to 10* iterations of a chaotic trajectory and the blank re- 
gion corresponds to the regular island. The straight lines 
represent three different families of period-three MUPOs (see 
text), and the different symbols correspond to specific MU- 
POs in one of these families, (b) Phase-space portrait of 
map ©-© for K = 2. We plot -0.5 < y < 0.5 in (a) and 
—0.5 < x,y < 0.5 in (b) for visualization convenience, (c) 
From bottom to top, the cumulative RTDs of the maps con- 
sidered in (a) and (b) (multiplied by a factor 10 for clarity). 
The upper curve is a straight line with slope 7 = 2. Inset: 
distance of a chaotic trajectory to the the border of the island 
in (a) during an event with recurrence time T = 335. 



regular island. The only relevant difference is that in this 
case there is no other family of MUPOs apart from the 
border of the island. 

In the stickiness of a chaotic trajectory, the trajectory 
approaches a family of MUPOs and follows a nearly pe- 
riodic dynamics for a long period of time before leaving 
the neighborhood of the MUPOs [inset of Fig.^c)]. In 
Fig. 12 we show the trajectories that stick to the MU- 
POs at the border of the regular island and remain in 
the neighborhood of the island after n' — 1, 2, 4 and 
1000 iterations of the continuous sawtooth map. In 
general, the closer to the island the longer it will take 
for the trajectory to leave. This stickiness is properly 
quantified in terms of the RTD for a recurrence region 
taken apart from islands. We have performed numeri- 
cal simulations for two different configurations present- 
ing sharply divided phase space: the continuous saw- 
tooth map with K = 3/2 [Fig. [Ha)] and map ©-(EI) 
with K = 2 [Fig. mb)]. As shown in Fig. HJc), in both 



cases the cumulative RTD is best approximated by a 
power law with scaling exponent 7 = 2. 

We have found similar results for other piecewise-linear 
area-preserving maps with polygonal islands, which were 
obtained from the maps considered above for different 
choices of the parameter K and from a different map 
studied in Ref. ^2§\. There are cases when the results 
described in this section and the theory of Sec. Ill C| do not 
apply, such as when the regular islands are elliptical and 
the outermost torus is quasi-periodic. Nevertheless, we 
observed that in many of these cases the exponent 7 = 2 
also fits the power-law tails of the numerically obtained 
RTD. 
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FIG. 2: Phase-space portrait of the continuous sawtooth map 
for = 3/2 showing the initial conditions of the trajectories 
that remain inside the dashed triangle for at least n' = 1, 
2, 4 and 1000 iterations of the map, respectively. The inner 
triangle corresponds to the regular island. We plot < a; < 
1.5 and —1 < x < 1 for visualization convenience. 



B. Mushroom billiards 

Billiards can be used as simple models in the study of 
Hamiltonian systems. Recently, Bunimovich introduced 
the so called mushroom billiards , which are billiards 
that have a single regular and a single chaotic ergodic 
region. A typical mushroom billiard is defined by a semi- 
circle (hat) placed on top of a rectangle (foot), as de- 
picted in Fig. |2ta) . The phase space is described by the 
normalized position x on the boundary of the billiard 
and angle 9 e [—0.5, 0.5] with respect to the normal vec- 
tor right after the specular refiection. The regular region 
corresponds to the orbits in the hat of the mushroom 
that never cross the dashed circle of radius r in Fig.|2Ia). 
The border between the regular and chaotic region of the 
mushroom billiard is therefore non-hierarchical, as shown 
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in figEIb). 

Mushroom billiards have two different classes of MU- 
POs, as illustrated in Fig.|3Ia). One of them corresponds 
to orbits bouncing between the parallel walls in the foot 
of the mushroom. Similar MUPOs are also found in many 
other billiards with parallel walls, such as the Sinai and 
stadium bilhards [lo, El ISHli • The other and more 
interesting class of MUPOs corresponds to periodic or- 
bits in the chaotic region that never leave the hat of the 
mushroom. In a previous study [l^ . we have shown that 
there is usually a complex distribution of these MUPOs 
close to the regular region and in contact with the chaotic 
component. The border of the regular island can be re- 
garded as MUPOs of this class. We have found similar 
MUPOs in other billiards with circular component, such 
as annular billiards j29i | . 

A relevant point concerning the stickiness in mushroom 
billiards is that the scaling exponent of the cumulative 
RTD is again 7 = 2, regardless of the control parame- 
ter < r/R < 1, as shown in Fig. ^c). In this case, 
the whole foot of the mushroom is taken as the recur- 
rence region in order to avoid the trivial parallel wall 
MUPOs. The injection and escape mechanism of the 
chaotic trajectories near the island are slightly different 
from those observed in the continuous sawtooth map be- 
cause in mushroom billiards there are escaping regions 
tangent to the island and the injection and escape oc- 
cur in a single iteration [inset of Fig. Qfc) vs. inset of 
Fig. EIc)] However, from a more fundamental point 
of view, the stickiness is remarkably similar in these sys- 
tems because in both cases the stickiness is mediated by 
MUPOs and the RTD has an exponent 7 = 2. Alto- 
gether, this suggests the possible existence of a universal 
scenario for the stickiness of chaotic trajectories in sys- 
tems with divided phase space, as considered in the next 
section. 



C. Scaling exponent: Theory 

We now derive the scaling exponent 7 = 2 for the cu- 
mulative RTD of two-dimensional systems with sharply 
divided phase space. Our theory apphes to the class of 
systems presenting one-dimensional families of MUPOs. 
This includes as particular cases some systems without 
regular islands, such as the stadium and Sinai billiards, 
whose cumulative RTDs are known to be governed by an 
exponent 7 = 2 [2^ l23l | . Most importantly, our results 
also apply to systems with mixed phase space, such as 
those considered in the previous sections. The theory re- 
mains valid when the orbit at border of the regular region 
is quasi-periodic and the first-escape region is tangent to 
the border. That is, if there are trajectories arbitrarily 
close to the island that move away in one or few time 
steps, such as in the mushroom billiards (18| . 

The essential features of the systems considered in 
the previous sections are captured by an area-preserving 
map M(x, 9) defined on the torus and that contains 




FIG. 3; (a) Typical musliroom billiard, defined by the geo- 
metric parameters (r, R, h). Two MUPOs are sliown in dotted 
lines, (b) Pliase-space representation of the semi-circular hat 
of the mushroom billiard with r/R = 0.5. (c) From bot- 
tom to top, RTDs for r/R = 0.6, 0.75, and 0.5 (multiplied 
by a factor 2 for clarity). The upper curve is a straight line 
with slope 7 = 2. Inset: distance of a chaotic trajectory 
to the border of the island during an event with recurrence 
time r = 300. 



a one-parameter family of MUPOs of period q. For 
concreteness, we assume that the family of MUPOs is 
{xi < X < X f, = 9o}. The phase space of map M(x, 9) 
is sketched in Fig.^a) and a possible configuration space 
in Fig. ^b). The following analysis does not depend on 
whether the MUPOs are in the chaotic sea or at the bor- 
der of a regular island. 

Consider small perturbations of a MUPO {xq, 9q): 

(i) If (a;', 9') = [xq + £xi 9o) and Xi < x' < Xf , another 
periodic orbit of the set of MUPOs is obtained. In 
this case M'^{x',9q) = (x',6'0), that shows that the 
perturbation neither grows nor shrinks. 

(ii) If {x',9') — {xo,9o + e), the perturbation in the 
9 direction does not grow. On the other hand, in 
the X direction the trajectory is not strictly peri- 
odic anymore and there is a displacement 6x every 
period q: M''{xo, 9') = {xq + 6x, 9'). 

Both effects (i) and (ii) have to be taken into account 
when a generic perturbation is considered: 

M'^ix', 9') = M«(xo + 6I0 + e) = {x' + 5x, 9'). (5) 

After q iterations, the same arguments used above 
for {x',9') apply to {x' + Sx,9'). We thus see that the 
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perturbed trajectory follows the dynamics remaining 
at a constant distance e from the family of MUPOs, until 
it travels Ax = x f — xq reaching the end x — Xf of the 
family of MUPOs (see Fig.^jl 30]. Wc note that Eq. © 
implies a linear growth of the perturbation in time, what 
is consistent with the marginal instability of the fixed 
point that forbids exponential growth of perturbations. 

The displacement 5x is related to the difference be- 
tween the frequency of the perturbed and unperturbed 
orbits, and can therefore be approximated linearly as 

5x = De, (6) 

in the limit of small e. For the continuous sawtooth map 
with K = 3/2, one obtains D = 6. In the case of bil- 
liards with parallel walls D = 21, where I is the distance 
between the walls. For MUPOs in circular like billiards, 
such as mushroom and annular billiards, D — 2qR, where 
R is the radius of the circle. 




FIG. 4: (a) Illustration of the dynamics of a perturbed 
MUPO {xo,0o + e) in the phase space (see text), (b) The 
corresponding dynamics in the configuration space of a bil- 
liard with parallel walls {q = 2). 



as a constant. From Eqs. Q and we obtain = 2 
for the power-law exponent of the distribution of escape 
times, or jtr — 1 for the cumulative distribution. This 
description is not valid when recurrences are calculated, 
because the initial conditions are chosen inside the re- 
currence region and thus away from the MUPOs. In this 
case, the convergence of p(e) to the invariant measure is 
much slower (algebraically for the stadium billiard 23]) 
and for any finite time p{e) — > for e ^ 0. However, we 
show in the Appendix that the scaling exponent for this 
second case can be derived from the first: the power-law 
exponent increases by -f-l when the initial conditions are 
taken away from the MUPOs, i.e., 7 = ■^tr + 1 (see also 
Refs. 0, 123, 1^ 113] ) ■ In systems with a one-parameter 
family of MUPOs, this leads to the asymptotic expo- 
nent 7 = 2 for the cumulative RTD, in agreement with 
our numerical results. 

Since every family of MUPOs contributes with the 
same exponent 7 = 2 asymptotically, the exponent does 
not depend on the possible presence of other families of 
MUPOs in addition to the one at the border of the reg- 
ular islands. Indeed, a large number of other families 
is observed in mushroom billiards |18| . a small number 
in the continuous sawtooth map and none in the case of 
map ©"O, and the scaling exponent is the same in each 
of these cases. 

It would be interesting to verify the generality of the 
exponent 7 = 2 33], for example, by investigating other 
systems presenting sharply divided phase space j34j . The 
stickiness through MUPOs described above resembles the 
mechanism underlying stickiness and anomalous diffusion 
in one-dimensional maps with marginally unstable fixed 
points 35]. Here we have considered two-dimensional 
systems and we believe that similar results hold true in 
higher dimensional Hamiltonian systems. 



The time a perturbed trajectory takes to reach Xf and 
escape from the dynamics (|SJ is given by 



for small e. In what follows, we see that this time is 
equivalent to the recurrence time if the initial condi- 
tions are chosen properly. Relation {Tj) shows that the 
smaller the perturbation the longer the time the trajec- 
tory takes to escape. The asymptotic distribution of es- 
cape times P{T) as a function of the distribution of per- 
turbations p{e) is given by 

The distribution p{e) depends on the choice of the initial 
conditions. 

For instance, choosing the initial conditions in the 
neighborhood of the family of MUPOs leads to a rapid 
convergence of p{e) to the invariant measure of the sys- 
tem. In this case, p{e) can be asymptotically regarded 



III. HIERARCHICAL PHASE SPACE 

Now we consider perturbations to systems with sharply 
divided phase space. This leads us to the problem of 
stickiness in Hamiltonian systems with the usual complex 
hierarchy of infinitely many KAM islands and Cantori. 
As a model system we consider the mushroom billiard 
perturbed by a magnetic field, a system that we refer to 
as the magnetic mushroom billiard and that allows for a 
direct comparison between the effects of hierarchical and 
non-hierarchical borders. 



A. Perturbation of non-hierarchical borders 

A common feature of the systems considered in the pre- 
vious sections is that their dynamics is piecewise smooth 
and presents abrupt changes. These abrupt changes, gen- 
erated by non-smooth functions / in map Q and sharp 
corners in the mushroom billiard, are responsible for the 
creation of sharply divided phase spaces. Generic per- 
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turbations of these systems are expected to smooth down 
the dynamics and introduce hierarchies of KAM islands 
and Cantori. Examples of such perturbations include to 
smoothen functions Q or Q in the case of piecewise- 
linear maps and soften the walls in the case of mushroom 
billiards. In the case of a billiard with charged particles, 
we can also perturb the system with a magnetic field, as 
studied below. 

Consider the mushroom billiard studied in Sec. Ill Bl 
subject to uniform transverse magnetic field B and con- 
sider the dynamics of charged particles within this bil- 
liard. Due to the Lorentz force, the charged particles 
move on circular orbits. We choose the charge of the 
particles and orientation of the magnetic field such that 
the trajectories are oriented counter-clockwise and have 
radius 

L (X 1 , (9) 

which is used as a control parameter. This parameter has 
to be compared with the geometric scales of the billiard 
defined in Fig. Ola) (in our simulations we use R — 2 and 
r = 1). The unperturbed mushroom billiard corresponds 
to L = oo. 
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FIG. 5: Magnification of the phase-space portrait of the mag- 
netic mushroom billiard at the border between the chaotic 
and regular regions for r/R = 0.5 and various values of the 
magnetic field. 

Previous works on magnetic billiards |361] have shown 
that the curvature of the trajectories often leads to the 
creation of KAM tori [stI Issj l in fully chaotic systems 
and chaotic regions |39l| in integrable systems. Mushroom 



billiards have both integrable and chaotic regions in the 
phase space and both effects are expected to take place. 
More interestingly, mushroom billiards also have MUPOs 
that are expected to undergo a transformation when the 
system is perturbed. Indeed, because the eigenvalues as- 
sociated to these orbits are real and have modulus 1, 
arbitrarily small perturbations are expected to generate 
elliptic or saddle points in the neighborhood of the reg- 
ular island of the unperturbed billiard. These effects of 
the magnetic field in the mushroom billiard are shown in 
Fig. 13 where a representative magnification of the phase 
space at the border of chaos is shown for different val- 
ues of L. The hierarchy of KAM islands and Cantori are 
clearly visible, providing evidence that the complete pic- 
ture of Hamiltonian chaos is obtained in magnetic mush- 
room billiards. 
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FIG. 6: Cumulative RTDs for the magnetic mushroom billiard 
with r/R = 0.5 and different values of the magnetic field. 
From bottom to top the lines represent: a power law with 
7 = 2, the numerical results for L = oo (shifted downward by 
two decades for clarity), L = 100 (shifted downward by one 
decade) and L — 50, and a power law with 7=1. 

The emergence of complex structures of KAM islands 
in the phase space influences the stickiness, as shown in 
Fig.|nifor RTD of magnetic mushroom billiards with L = 
100 and L — 50. Comparing these distributions with 
those of the unperturbed system (L = oo), we note the 
presence of fluctuations around a slower power-law ten- 
dency (7 < 2). This result indicates that, as intuitively 
expected, a hierarchical border sticks the trajectories in a 
more effective way than a non-hierarchical border. While 
the presence of a single family of MUPOs in a hierarchi- 
cal phase space would be enough to guarantee 7 < 2, 
we observe that, generically, all the families of MUPOs 
disappear. One could expect that the outermost torus 
of a regular island, which is marginally unstable, could 
play the role of the MUPOs described in Sec. ^1 How- 
ever, there is usually an infinite number of Cantori that 
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accumulate near the island invalidating relations jSJ and 
© and thus the derivation of the exponent 7 = 2. In 
the next section we study carefully the effect of such a 
hierarchical border on the stickiness. 



B. Hierarchical phase-space scenarios 

We now investigate the origin of the oscillations in the 
RTDs shown in Fig. El 

We focus initially on the parameter L = 50. For this 
parameter, many KAM tori are destroyed but the chain 
of islands and Cantori are still clearly visible in the phase 
space, as shown in Fig.[7| The different density of points 
seen in Fig. Cfa) is related to the presence of chains of 
islands and Cantori acting as partial barriers |^ to the 
transport in the 9 direction. In order to associate the 
presence of these barriers to the RTD, we study the min- 
imum distance between the trajectory and the main is- 
land before the trajectory leaves the neighborhood of the 
island and visits the recurrence region. In our simula- 
tions we use the minimum collision angle of the tra- 
jectory as a measure of the distance because the barriers 
mimic the original tori and have approximately constant 
6, and we take the foot of the mushroom billiard as the 
recurrence region. The fraction of events that have a 
minimum angle 9 is defined as g{9)d9 = rio/rj, where rjg 
is the number of recurrences that have a minimum angle 
in the interval [9, 9 + d9] and rj is the total number of 
recurrences. Numerical results for g(0) with L = 50 are 
shown in Fig. [7fb). The function g{9) goes to zero at 
the angles that correspond to the position of the barri- 
ers because the trajectories that manage to pass a bar- 
rier quickly spread throughout the next chaotic region. 
From the behavior of g{9) in Fig.[7{b), we can resolve 5 
different regions limited by these barriers. To associate 
these regions with the RTD, we label all the recurrence 
events from (1) to (5) according to the number of re- 
gions the trajectory penetrates before returning to the 
recurrence region. The RTD of each of these groups of 
recurrence events are shown in Fig.|3c). The RTD of all 
the events corresponds to the sum of these partial RTDs 
and is shown in the same figure (upper solid curve). The 
partial RTD of each region (l)-(5) presents a relatively 
peaked maximum followed by an exponential decay. Ac- 
cordingly, most of the orbits that have the same recur- 
rence time T penetrate the same number of barriers [note 
the logarithmic scale in Fig. [7{c)] . These results indicate 
that, for T < 10^, the stickiness is dominated by the pri- 
mary chain of barriers around the main regular island, 
that is, the contribution of barriers associated to sec- 
ondary islands is negligible. These results also show that 
the oscillations observed in the RTD around the power- 
law behavior are intrinsically associated to the presence 
of the barriers in the phase space. 

These stickiness properties agree well with the predic- 
tions of the model proposed by Hotter et al. in Ref. . 
In that paper, the hierarchy of Cantori is modeled by a 



(a) (S) 




T 



FIG. 7: Analysis of the magnetic mushroom billiard with L = 
50: (a) phase-space magnification at the border of chaos; (b) 
fraction g[6) of recurrences that have 6 as their minimal angle; 
(c) the RTD of all trajectories (upper solid curve) and the 
RTDs of the trajectories in regions (l)-(5) of (b) (lower solid 
curves). The lower curves in (c) are divided by 10 for clarity. 



chain of coupled hyperbolic systems, where each hyper- 
bolic system models the area of the phase space limited 
by successive Cantori. One of the strengths of this model 
is that it predicts not only the asymptotic behavior of the 
non-hyperbolic dynamics around KAM islands but also 
the finite-time dynamics assessable in numerical simula- 
tions and experiments. The model predicts that the sur- 
vival probability of particles in the neighborhood of KAM 
islands fluctuates around a power law and is composed 
of a sum of exponentials associated to the Cantori. Our 
results in Fig. [7| show that this behavior is indeed present 
in real Hamiltonian systems. As shown below, this pic- 
ture changes when secondary structures of the hierarchy 
are relevant. This more general stickiness scenario is ob- 
served in the mushroom billiard for larger values of the 
magnetic field (e.g., L — 10). 

In Fig. |S1 we show the same as in Fig. [7| for the pa- 
rameter L = 10. The effect of the primary barriers is 
still important, as shown in Fig. Elb) where these bar- 
riers correspond to zeros of g{9). However, as shown 
in Fig. |S{c), the partial RTDs corresponding to regions 
(l)-(5) exhibit a power-law rather than an exponential 
decay. For instance, the RTD of trajectories belonging 
to region (2) exhibits an approximate power-law decay 
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FIG. 8: (Color Online) Analysis of the magnetic mushroom 
billiard with L = 10. (a) Phase-space magnification with two 
typical sticking trajectories with recurrence time T ~ 8 10'': 
trajectory 1 sticks near the upper island and trajectory 2 fills 
the chaotic region, (b) Fraction g{9) of recurrences that have 6 
as their minimal angle, (c) The RTD of all trajectories (upper 
solid curve) and the RTDs of the trajectories in regions (1)- 
(5) of (b) (see legend). The lower curves in (c) are divided by 
10 for clarity. 



that makes these recurrence events dominant not only for 
small times (10 < T < 500) but also for very large times 
(T « 10^). On the other hand, the RTD of events associ- 
ated to region (4) does not dominate the (total) RTD at 
any time. The slower decay of the RTD of region (2) is 
a consequence of the stickiness to the chain of secondary 
islands shown at the top of Fig. IS^a). In this figure, 
we show two representative trajectories with recurrence 
time T « 8 10*. The first (trajectory 1) penetrates only 
two regions and sticks to a secondary island. The second 
(trajectory 2) penetrates five regions and approaches the 
main island. In the context of stochastic models p3.l40l|. 
asymptotic effects of secondary islands can be accounted 
for by the Markov-tree models 0. In a deterministic 
framework, the full hierarchy of islands can be accommo- 
dated within a chain model of non-hyperbolic systems or 
a tree model of hyperbolic systems, which are straightfor- 
ward generalizations of the model introduced by Motter 
et al. 



IV. CONCLUSIONS 



We have studied the stickiness of chaotic trajectories in 
Hamiltonian systems with sharply divided phase space, 
which are characterized by non-hierarchical borders be- 
tween the regions of chaotic and regular motion. The 
stickiness occurs through the approach to one-parameter 
families of MUPOs in contact with the chaotic region. 
The main characteristics of this stickiness scenario are 
the exponent 7 = 2 for the power- law decay of the cumu- 
lative RTD and the long intervals of regular motion at a 
constant distance from families of MUPOs. Dynamical 
systems described by this scenario include mushroom bil- 
liards and various piecewise-linear area-preserving maps. 

Generic perturbations applied to systems with sharply 
divided phase space destroy the MUPOs and introduce 
hierarchies of regular islands and Cantori. We believe 
that these perturbations can serve as a new paradigm to 
the study of stickiness in generic Hamiltonian systems. 
Using as an example mushroom billiards perturbed by 
a transverse magnetic field, we characterize two differ- 
ent scenarios of stickiness in the presence of perturba- 
tions. For small perturbations, the stickiness is domi- 
nated by the primary chain of Cantori, which work as 
partial barriers to the transport around the main regular 
island. In this case, the RTD is composed of a sum of 
exponential distributions associated to the probability of 
crossing each of these barriers. For increasing perturba- 
tions, the primary barriers weaken while the secondary 
islands and the corresponding sticking regions grow. For 
large perturbations, the stickiness of the secondary struc- 
tures becomes relevant and the exponential components 
of the RTD are converted themselves into power-law dis- 
tributions. This provides direct evidence of the effects of 
Cantori structures at finite times, in strong support of 
the model introduced in Ref. and its generalizations. 

The asymptotic behavior of the RTD, which has been a 
matter of considerable recent debate 5,6], cannot be re- 
solved alone by numerical experiments. Our simulations 
suggest that the hierarchical structures enhance the stick- 
iness of non-hierarchical borders, what would lead to an 
upper bound 2 for the scaling exponent 7. This upper 
bound is guaranteed when the phase space has one or 
more families of the MUPOs described in Sec. ^ How- 
ever, in general, this numerical evidence of upper bound 
should be taken with caution because one cannot neglect 
the possibility that the hierarchical structures will reduce 
the stickiness for asymptotically large times. For general 
Hamiltonian systems, even the question of whether the 
oscillations in the RTD vanish asymptotically, giving rise 
to a well defined power-law exponent, is a problem yet to 
be settled. Our results provide an answer to this ques- 
tion for an important class of Hamiltonian systems with 
sharply divided phase space. 
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Appendix 

Consider a uniform distribution of initial conditions in 
the neighborhood of a sticking region of the phase space, 
as usually studied in problems of transient chaos, and 
consider the time it takes for the corresponding trajec- 
tories to escape to a pre-defined region away from the 
sticking region. The distribution S{t) of escape times 
longer than r is proportional to the measure ^(r) of the 
region of the phase space to which the trajectories stick 
for a time longer than r. Due to the ergodicity, 

S{t) (X M(r) - ^, (10) 

where It is the total time spent inside the sticking region 
and t is the total observation time. 

On the other hand, in the study of recurrence prob- 
lems, one usually initializes a single trajectory in a re- 
currence region away from the sticking region and com- 



putes the time T the trajectory takes to return to that 
region. If the trajectory is followed for a long time t, the 
cumulative RTD is 

Qir) = ^, (11) 

where A'',- is the number of recurrences with time T > t 
and A'^ is the total number of recurrences observed in 
time t. The relation between the times in Eq. IjlOII and 
the number of recurrences in Eq. (|11|) is given by j5|] 

t ^N{T), (12) 

- Nr T, (13) 

where (T) is the average recurrence time. Altogether, 
this leads to 

Q{r) - ^ . (14) 

T 

In particular, if the escape times follow a power-law dis- 
tribution S(t) ^ r then the cumulative RTD is 
Q{t) T where 

7 = 7«r- + 1 • (15) 

An equivalent relation was obtained in Ref. |32l | for 
chaotic scattering problems, another example where the 
trajectories are initialized away from the sticking region. 
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